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LINEAR ACCELERATION GUIDANCE FOR LUNAR 
LANDING AND LAUNCH TRAJECTORIES 

SUMMARY 


An approximate analytical solution to the problem of maneuvering a 
spacecraft to reach specified end conditions by finite thrusting in the 
vicinity of the moon is presented. The analysis includes all three 
dimensions and the acceleration terms in the equations are linearized. 

The examples presented are limited to the planar case. Although in 
general the thrusting scheme is variable, a constraining relation is 
introduced to allow burning time for constant thrust to be calculateu. 
Similarly, a constraining relation may be introduced which allows a burn- 
ing time for a constant thrust angle to be calculated. 

The solution is versatile in that it allows variable thrust, constant 
thrust, and constant thrust angle trajectories between specified end con- 
ditions. 


INTRODUCTION 


In order that a spacecraft such as the Apollo Lunar Excursion Module 
(LEM) can be guided during landing, launch, vJ>ort or rendezvous to c. set 
of specified end conditions in a manner that is near fuel optimal, a set 
of guidance equations must be mechanized on board the spacecraft that 
will predict the necessary thrust and/or thrust angles. This paper will 
present the derivation of a set of equations that are suitable for this 
task from the standpoint of guidance and fuel optimum performance. 

In obtaining a solution for the guidance equations, it is desirable 
that several criteria be met: 

(1) The equations must be computationally simple. 

(2) The equations must be suitable for use in as many operational 
modes as possible. 

(3) The equations must yield a solution the;t is near the fuel 
optimum. 

In order to obtain such a set of guidance equations, tne first step 
was choosing a suitable approximation to linearize the equations of 
motion. To accomplish this, it was assumed that the change in altitude 
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of the spacecraft is small compared, to the initial radius of the space- 
craft. This is similar to assuming a constant gravitational field. The 
next step was to soxve the two- point boundary value problem explicitly. 

The approach taken in this paper in solving the problem is to pre- 
scribe that the radial, tangential, and out-of-plane components of the 
acceleration vary linearly with time. The equations of motion can then 
be solved in closed form. 

Six parameters are introduced in the three acceleration components 
which can be determined in closed form in terms of the six specified end 
conditions. These parameters constitute the guidance equations which 
will always insure that the spacecraft's trajectory will meet the speci- 
fied terminal conditions. 

For this acceleration scheme, the equation for the thrust angle is 
the same as that for the fuel optimum for a flat central body or constant 
gravity field approximation where the final position and velocity are 
constrained (ref. l) . 

It is also shown that by introducing other constraining relations, 
constant thrust and constant pitch angle trajectories can be generated. 


LIST OF SYMBOLS 


A^> Ag, Ay A^, A^ 


constants defined on page 


a 

b 


c 



F(t) 


condition for constant thrust 


G(t) 


condition for constant pitch angle 


g 


E 


52.2 yy-, constant relating mass to weight in 

sec 

earth pounds 



altitude, feet 


t 
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•V:' 


\ : 
-X. 


>V 
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■r - 
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I 

;p 



Kj, 


m 

R 

r, 0 , ^ 

T 

t 

v 

V 

c 

W 

X 

z 


specific impulse, seconds 

arbitrary constants of integration related to 
initial conditions 

mass 

radius of attracting body, feet 

polar coordinates defined in Sketch 1 

thrust, pounds 

time, seconds 

velocity, feet per second 

characteristic velocity, V , feet per second 

weight of spacecraft, pounds 

arc length ox range along surface of 
attracting body, feet 

out-of-plane range, feet 


a initial level of applied acceleration in a 

2 

given direction, ft/sec 

0 thrust azimuth angle, degrees 

0 (with subscripts) rate of change of applied acceleration in a 

given direction, ft/sec^ 

S( ) denotes variation 

G angle between thrust vector and local 

horizontal; the pitch angle, degrees 


* 


i ■ 
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V 


(') 

( )' 


Subscripts 

0 

1 , 2 , 3 , 4 , 5 
r 

T 

0 

♦ 


universal gravitational constant times mass of 

5 2 

attracting body, ft '/sec 

T/g„I m . seconds" 1 
JCj sp c 

burning time, seconds 

derivative of ( ) with respect to t 

derivative of ( ) with respect to t 

conditions when t = 0 
denotes different A's and K's 
pertaining to radial direction 
conditions when t = t 

pertaining to ci: cumferential, or x, direction 

pertaining to direction perpendicular to r, 0 
plane 

DERIVATION OF EQJATIONS 


The equations of notion of a thrusting spacecraft in a gravitational 
field are: 

r - r(0 2 + ijr 2 sin 2 0) + ~ 

r 


d , 2. . 2*2 A 0 
rr (r 0) - r f sin 0 cos 0 = r -£• 
at m 


d / 2* \ 2 * * i|/ 

tt- ( r V sin 0) + r i|r 0 cos 0 = r — 

GLt XU 


(1) 

( 2 ) 
( 3 ) 


See reference 2. 



These equations may be solved analytically by prescribing linear 
acceleration components: (See sketch 1.) 

T 

— = — sin 0 = a + 3 t (M 

mm r r 

T 

— = — cos 0 cos 0 = a + 0 t ( r ) ) 

mm 0 0 

T 

— cos 9 sin 0 ~ a + 0 . t (l) 

mm ^ t 

and assuming thau the change in altitude ;r» small compared to the initial 

radius r , and that the. out-O'-; -ant angle and its time derivative, 
o 

and are small. 

In equation (2), assume that r = i = constant that: 

2 * ? d 2 * 

r ijr~ s^n 0 cos 0 < < ™ (r 0 ). This equation, using •" > then becom< 


U 

dt 



r (a 
o ' 0 


+ V° 


Integrating this equu ; .. . on twice, 


0 * 


0 - 



, * 0 t + » ; 


+ K, t + x a t" r 
1 ^0 



( 7 ) 

(8) 


Before applying initial conditions to ( 7 ) and (8), make a change 
of variables by defining the range X - R0. Equations (7) and (8) then 
become: 
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where and are constants of integration. 

•2 2 .2 

Assuming that r = r = constant and that ^ sin 0 < < 0 , 

o 

equation (1) becomes, 

r - r 0^ + = a + 8 t 

o 2 r r 


( 11 ) 


# 2 

Square both sides of equation (7)., and then substitute r 0 from 

on, ° 


r 

o 


- a - p t 
2 r r r 


j 


-2 k 2 * + y 2 + * 2 

n L 

H 2 * 4 ] 


+ apt + 
0 0 


Solving for r, 

r 


0 ' / 


r 4r 


( 12 ) 


Now, r = h 4 R, where h is the altitude of Jie spacecraft and R 
is the radius of the attracting body. Using r = "ft and by defining, 

“i 2 

* r- ♦ «r - -“5 

o r 

o 


2K,a 

Ag - -j- 2 - * e r 
o 

Aj - i- (a, 2 ♦ V 0 ) 


I 


A, 


a B 
« 0 0 



equation (12) becomes: 


f ‘ =A l 


2 4 

+ Agt + A^t + A^t' + A^t 


(15) 


Now integrate (13), 



J 


+ lv 2 


+ F A 3 t? + 4 A 4 t 


+ i A.t" 
5 ? 


(14) 


12'' 3 1 4 

h = K 4 + K^t + ^ k^t + g Agt + 12 V 


+ A V 5 + 30 V 6 


(15) 


where K and K. are constants of integration. 

3 4 

Define the out-of-plane range, Z = r sin 0 sin «jr 
For small angles f ? and assuming r = r Q = constant this becomes. 


Z = ♦ r Q sin 0 


and the out-of-plane velocity becomes, 

- Z » ♦ r Q sin 0 + * 0 r Q cos 0 

Substituting this last expression into (3), and again assuming 

r = r = constant, (3) becomes, 
c 

A. (Z - r *gl cos 0) + r Jr0 cos = cl, + P..t 
dt o ^ i V 



Differentiating the first term and using the definition X = R0, the 
last equation reduces to: 

2 + T' i '« (iccos0) ‘V V (l6i 

Since ^ is a small angle and if X does not change too rapidly 
with time, the second term on the left may be neglected and (l6) becomes 


Z = a, + 8 t 

+ t 

Equation (17) may now be integrated twice to give: 


i ‘ «5 + V + 5 V 2 


7. = K 6 + Y + i a/ + \ Y 3 


(17) 


(18) 

(19) 


Equations (9), ( 10 ), (l 4 ), (lp)> ( 18 ), and (19^ are the equations 

of motion of the spacecraft. Now apply the initial conditions to 

evaluate K, through K,. When t= 0 , X=X,X = X,h = h,h = h, 

1 # 6 o o' o' o 

Z = Z q , Z = Z q . The constants of integration are seen to be: 
r r 

KL. = X , K. = —• , K, = h , K, = h , K,. = Z , and K. = Z . 

1 R o’ 2 R ’ 3 o’ k o’ "5 o’ 6 o 


The equations of motion are now determined except for the thrust 
parameters ct^, 0^, a^, 0^, a^, and 0^, which determine the thrust 

history; that is, the thrust magnitude and direction. The six equations 

of motion may be solved for these thrust parameters in^ terms of final 

conditions. When t=T, X-X,X=X, h = h, h=h, Z = Z , and 

T T T T T 

Z = Z^ . After considerable manipulation the equations for the thrust 
parameters become: 


O r -1 

“0 (X o-V + 2 *o + *tJ < 20 > 

p 0 • 7 r[f < x o - V + K * * T ] < 21 > 


m 


o i A_ A. 

a = - -§• X + — ( *h - h ) - -^■T'-r^T - -p- -r' ? 

2 _2 o T T o 2 3 4 


r o R 




( 22 ) 


o 12 /, , » 6 ,♦ • v 

P = - <zo. — -r —r \n - n / + -s in + n ; - a, t 

r 0 R 3 o r 2 ' o t 3 

T T 


9 2 4 3 

A, T - ~ A, T - — A 

3 10 4 5 5 


(23) 


2 

a = - - 

i i t 


[1 (Z - Z ) + 2Z + Z 1 
It o t o t 

13 = 4-T “ ( z - Z ) + Z + zl 

♦ t 2 |_t c t 7 o tJ 


(24) 


(25) 


For a specified set of initial and terminal conditions the boundary 
value problem is now completely solved. 


The thrust azimuth angle, j3, is given at any time by dividing 
equation (6) by (5), 

T, An a + p t 
tan 3 - JL » JL J- 

V“ a 0 + v 


( 26 ) 


The thrust pitch angle, 0, is given at any time by dividing (4) by 
(5) and using (26), 


tan 0 = 


V 


a + 3 t 


cos 3 T 0 /* cos 3 a 0 + 3 0 t 


(27) 


T 

The total acceleration, — , is given at any time by, 


/ V2 /T \ 2 /T \ 2 /tA 2 

It) . i-y Ji) 

\ m J 1 m } V m / \ m / 


= (» r + 0 r t) 2 + (a 0 + 3 0 t) 2 + (a + + 3^,t) 5 


« 
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I = + bt + Ct 2 (28) 

where, 

2 2 2 
a = a + a + a, 
r 0 + 

b » 2( ar P r + a 0 3 0 + a^) 

2 „ 2 „ 2 
c = 3 + + P 

K r 0 y 

The total acceleration in (28) is now used directly in the definition 
for characteristic velocity, or performance index, of a given trajectory. 


v c ■ l (l) at - / y* + bt * <* 2 ~ « 


(29) 


Performing the integration equation (29) becomes: 


u ^ ac ~ o n 
c 8 c 5/2 


”^a + br + ct 2 + ^ ^~r 


. 1 b 

v a + 2 7T 


+ ?£ ic + ' b -yj* + bT + ct2 - (30) 

Since a, b, and c are dependent only upon the initial and final 
conditions and the burning time t, the characteristic velocity is related 
only to the burning time once the end conditions are specified. In most 
cases, it is desirable to minimize subject to certain constraints. 


THRUSTING MODES 


Giver, a complete set of initial and terminal conditions, there 
exists an infinity of trajectories that will satisfy the end conditions. 

If a burning time is specified, the trajectory is uniquely defined. In 
some cases, a burning time will be chosen that will minimize V . Two 

Important aite. n»te cases of course are, i l) choosing the burning time 
so that a constant thrust magnitude trajectory can be obtained, and 
(2) choosing the burning time so that a trajectory with a constant value 
of 0 is obtained. 

Case a: Minimum V . To obtain the burning time for minimum V , 

c c 

the derivative of V with respect to burning time may be set equal to 

c 

zero, and the resulting relation solved for t. This value of t must 
then be further investigated to determine if it yields a maximum, minimum, 
or an inflection point. This could be done but would be difficult, since 
V c is a very complicated function of t. It is much easier from the 

standpoint of mathematical complexity to solve (29) numerically by just 

evaluating it for a wide range of t and then looking for the minimum 

on a plot of V c versus t. Once the optimum value o: t is obtained, 

a , 3 , a , 3 , a,, 0 . can be calculated from (20) through (2p); 
r r 0 0 V \|r 

then the trajectory with minimum V c based on these value- is calcu- 
lated from ( 9 ), (10), ( 14 ), (15), (18), and (IS) • 

Case b: Constant Thrust Magnitude . Equation ( 28 ) may be shown to 

yield a constant thrust magnitude by comparing it to the acceleration 
relation for constant thrust. The acceleration for constant thrust may 
be obtained by considering the basic relation between thrust and rate 
change of mass. 


For constant thrust, 


Equating this to ( 28 ), 


T = " g E h V A = ™ 


/, t x • T T / m o 
m = m Q (1 - v t), v = - = I — 


V - 


. T Ac 


g E I sp 


/ — 
- « Ma + bt + ct 
m V 



(3D 




r 
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For a specified set of end conditions equation ( 31) may be solved 
for the burning time t, since a, b, and c are functions of t, at any 

time t. In particular at time t = 0, (31) becomes ^ a = T/m . 


2 2 2 

But since a = a + a + a, , the relation to be solved for t is, 

(A r '‘Ir ' 7 


a = 


2 2 2 ,2 • / t \ 2 2 

a 0 +a r +a * - ^V m o> ‘ (vTj ®E 


(32) 


To determine the proper constant thrust level, equation (32) may be 

solved simultaneously with the characteristic velocity equation for 

constant thrust, V » - g_ I In (l - vt)* 
c a sp 


Equation (32) must be solved for t by am iterative process. 
Define the function, 

V2 „ 




8x 


(33) 


and its derivative. 


£ F (*> - F '< T ' - + a r a P 


(34! 


The iteration process, 

r n+J 


f (t„) 


T ' * T n‘F7?P 


(35) 


is then set up to determine r. 

Equation (35) has been observed to converge quite rapidly for the 
constant thrust burning time. 

Case c; Constant 0 . Assuming 3 is small so that cos 0 « 1, 
the time derivative of equation (27) for 6 is. 


§ » j £ .2L — r P cos 0 

For 0 to be constant, 0 » 0, 60 (36) yields, 


156) 


“A - V* 


(37) 


\}4*l 



*%?• * 



i 

> 

j 


I I 


■Jr ’ 1 ' 1 


Hi 





1 ; 


How define the function, 


G(t) 


a,P_ 

0 r 


a B 
r 0 


= 0 


and its derivative, 

G'(t) = (a, p' + a'p ) - (a p' + a'8 ) 
v v 0 K r 0 r' v r 0 ff' 

The iterative solution for the burning time for constant 9 is, 

0(T n ) 

T n+1 ‘ T r. " OTT) (53) 


USE OF THE ANALYTICALLY DERIVED THRUST PARAMETERS FOR GUIDANCE 


Equations (20) through (25) give analytical expressions for the 
thrust parameters. They were derived on the basis of a model that 
approximated the inverse square gravity field. These equations may be 
considered as guidance equations which will guide the spacecraft during 
various maneuvers . The performance of these guidance equations can be 
evaluated by a simulation that is programed on a digital computer. A 
trajectory program which numerically solves equations (l), (2), and 
(5) to predict the position and velocity of the spacecraft at any time 
is a basic component of the simulation. A guidance program which contin- 
ually computes the thrust and its direction by solving equations (20) 
through (25) is the other basic component of the simulation. The 
instantaneous thrust vector computed by the guidance program depends 
only upon the present position and velocity of the spacecraft and the 
specified terminal, conditions . 


The problem of combining the two programs for variable thrust is 
straightforward. Both secs of equations are given the same initial 
conditions. The final conditions are also substituted into the thrust 
parameter relations, equations (20) through (25). p^, a , p , a 


y 


P ( are computed with these conditions and are used to calculate 9, p, 

and (T/m) from equations (26), (27), and (28), for the first step. The 
output, x, x, h, h, Z, 2, of the first integration step are then used 
as the initial conditions in equations (20) through (25) to compute 9, 
P, and (T/m), for the next integration step. The calculation thus pro- 
ceeds to the final conditions of the problem. 


If constant thrust or constant 8 is desired the burning time is 
updated by equations (55) or (58), just as 8, p, and (T/m) are in the 
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variable thrust case discussed above, t is updated prior to each 
integration step using the output from the previous step as the initial 
conditions. The computation procedure can best be visualized by re 
ferring to the block diagram in Sketch 2. 


EXAMPLES OF CALCULATIONS 


Several examples will now be given to illustrate the various modes 
in which the computation scheme may be operated. In all of the following 
examples, the calculations were limited to the two-dimensional case. 

(a) Variable thrust with minimum characteristic velocity: 

This mode is illustrated by a descent trajectory from 
approximately 50,000 feet to 10,000 feet above the lunar surface. 

The burning time was 274 seconds and was calculated by solving 
d V 

-~2. numerically. The V n was 4,833 feet per second and the total 

range was specified as 169. 7 n.mi. Time histories of the variables 
are shown in figures (la) through (if). 

(b) Constant thrust : 

This mode is also illustrated by a lunar descent trajectory^ 

The initial and final conditions are the same as in the variable 

mode discussed in (a) above. The burning time was 299*5 seconds 

and the V was 4,883 feet per second. The time histories of the 
c 

variables are shown in figures (2a) through (2e) . Also plotted for 

comparison is a trajectory having the same boundary conditions, 

that was calculated from a calculus of variations optimum program. 

This optimum trajectory had a burning time of 300. 90 seconds with 

a V of 4,883 feet per second. The V of the guided constant 
c c 

thrust trajectory was less than that for the optimum because of 

the slight difference in burning time between the two. It should 

be noted here that the two trajectories flown are different, except 

for the boundary conditions. The thrust-to- initial-weight ratio 

was 0.4. 

(c) Constant thrust launch with an error analysis for variations in 
the thrust vector: 

This trajectory had an initial altitude of 1,000 feet and an 
initial vertical velocity of 100 feet per second. The final 
conditions were pericynthion altitude of 50,000 feet, with orbital 
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velocity of about 5>55S feet per second. The total range was 92.92 n.mi. 
The guided trajectory is compared with the optimum in table I. Additional 
calculations of this trajectory were made perturbing independently the 
thrust magnitude T and thrust direction 0 by small amounts . The 
effect of these perturbances are given in table II. It is evident from 
this table that the end conditions are not seriously affected by 1.0 to 
5.0° misalinements in thrust direction, and 1.0 or 2.0 percent variations 
in thrust magnitude. The end condition most affected was the range. 


CONCLUDING REMARKS 


An approximate analytical solution to the problem of maneuvering a 
spacecraft to reach specified end conditions by finite thrusting has 
been presented. The feasibility of using this solution to guide the 
spacecraft to a specified final state has been demonstrated. Further, 
it has been shown that by the introduction of certain constraining 
relations that constant thrust and constant thrust angle trajectories 
may be flown. The guidance equations obtained from the analytical 
solution have been shown to be computationally simple; versatile, in 
that several operational modes arc possible; and economical in that they 
guide the spacecraft along a trajectory that is near the fuel optimum. 
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TABLE I.- COMPARISON OF FINAL CONDITIONS BETWEEN AN 

OPTIMUM CONSTANT THRUST LAUNCH AND A GUIDED 

LAUNCH; T/w * 0.6; I = 315 sec. 

* ' o sp 


Initial Final Conditions 

Conditions Optimum Guided 


x, feet 

0 

564,573.4 

562,982.0 

h, feet 

1,000 

49,957.94 

H~>, 956.62 

-, fps 

0 

5,558.02 

5,558.93 

h, fps 

100 

0.08257 

0.26230 

v f * s I 


5,775.25 

5,778.12 

1 

T f sec 

i 


227.916 

228.00 









EFFECT OF ERRORS IN PITCH ANGLE AND THRUS 


























Fig H: continued 

(e) Thrust angle as a function of time 






(a) Range as a function of time 
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